Measurement back-action on the quantum spin-mixing dynamics of a spin-1 

Bose-Einstein condensate 
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We consider a small F = 1 spinor condensate inside an optical cavity driven by an optical probe 
field, and subject the output of the probe to a homodyne detection, with the goal of investigating 
the effect of measurement back-action on the spin dynamics of the condensate. Using the stochastic 
master equation approach, we show that the effect of back-action is sensitive to not only the mea- 
surement strength but also the quantum fluctuation of the spinor condensate. The same method is 
also used to estimate the atom numbers below which the effect of back-action becomes so prominent 
that extracting spin dynamics from this cavity-based detection scheme is no longer practical. 

PACS numbers: 03.75.Mii, 03.75.Kk, 03.65.Ta, 42.50.Pq 



I. INTRODUCTION 



The ability of optical dipole traps to simultaneously 
cool and trap ground-state atoms in different magnetic 
sublevels paved the way for the experimental realization 
of spinor Bose-Einstein condensates (BEC) [l[ where the 
liberation of spin degrees of freedom has added new and 
exciting possibilities in the study of BEC Q • Of particu- 
lar relevance to the present work has been the extensive 
study of s pill dynamics both in theory and in exper- 
iments (7h11| , which aims to understand how condensate 
populations exchange coherently among different internal 
spin states as well as to explore the potential of these spin 
oscillations as probes to the intriguing physics underly- 
ing spin-dependent collisions. So far, such studies have 
been carried out mainly in spinor condensates with suf- 
ficiently large numbers of atoms where the measurement 
based on the standard absorption-imaging technique has 
established that the spin dynamics agrees well with the 
result predicted within the framework of mean-field 
theory 

The use of the spinor condensates with relatively large 
atom numbers, unfortunately, renders it virtually im- 
possible to observe bcyond-mean-field quantum effects, 
which, besides being fascinating by their own rights, are 
thought to be responsible for exotic physics in highly cor- 
related systems. In small spinor BECs, quantum fluctu- 
ations of atoms are expected to play a more prominent 
role. Hence, their spin dynamics should be governed by 
spin-mixing of quantum mechanical nature, which is re- 
sponsible, for example, for the generation of squeezed 
collective spin states and entangled states [l^- However, 
the absorption-imaging technique, which works well with 
large condensates, is no longer an effective detection tool 
for condensates with small numbers of atoms as a result 
of the much reduced signals. Recently, owing to the ex- 
perimental realization of strong coupling between ultra- 
cold atomic gases and cavity [l3j , spinor BEC confined in 
a sin gle- mode cavity has received much theoretical atten- 
tion |14| - [l7| . In particular, cavity transmission spectra 



have been suggested as candidates for probing the quan- 
tum ground state [l^ or the quantum spin dynamics [TBj 
of a spinor BEC. 

In all these proposals, one aims to learn the condensate 
dynamics from the photons landing on the photoelectric 
detector. The detection is a process, according to the 
Copenhagen interpretation of quantum mechanics, that 
projects the system (condensate -|- photon) to one of the 
states that are eigenstates of the observable being mea- 
sured. This causes the condensate dynamics to be dis- 
rupted in a random fashion and is the underlying physical 
mechanism behind the measurement back-action. Moti- 
vated by the possibility that such a back-action may be 
significant for a small condensate, we investigate, in the 
present work, the effect of the back-action on the quan- 
tum spin-mixing dynamics of a small spinor BEC subject 
to a homodyne detection as shown in Fig. [T] (sec next sec- 
tion for a detailed description). We take the stochastic 
master equation approach which emulates the experimen- 
tal process where many runs of continuous measurements 
must be performed before one can arrive at the quantum 
mechanical average of a dynamical observable. We show 
that the effect of back-action is sensitive to whether the 
condensate is in the ferromagnetic or in the antiferro- 
magnetic ground state. We point out, in connection with 
the recent interest in the dynamics of small spinor con- 
densates, that there is a limitation to this cavity-based 
detection scheme: for sufficiently small condensates, the 
number of experimental runs required to faithfully ex- 
tract the internal spin dynamics can be unrealistically 
large; we estimate the atom numbers above which this 
cavity-based homodyne detection scheme is experimen- 
tally feasible. 

The paper is organized as follows. In Sec. II we present 
a measurement model and derive the stochastic master 
equation (SME) describing the evolution of the spinor 
BEC. In Sec. HI, we first analyze the measurement back- 
action in two-atom case where some typical effects, such 
as thequantum Zeno effect (QZE) and quantum diffu- 
sion [3 [13: are clearly shown, and then extend the 
analyses to A'^-atom cases by numerical simulation with 
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realistic experimental parameters. We show different re- 
sponses of ferromagnetic and antiferromagnetic ground 
state to measurements, and measurement-induced deco- 
herence effects in quantum spin-mixing dynamics. In 
Sec. IV, the measurement outcomes are given by aver- 
aging over the photodetector currents of repeated mea- 
surements. Finally, Sec. V concludes the paper. 



II. MODEL AND HOMODYNE DETECTION 
SCHEME 

Figure [1] is a schematic of our model made up of three 
parts: a F = 1 spinor BEC, a driven ring cavity, and 
a homodyne detection arrangement. The BEC is as- 
sumed to be sufficiently small (with less than 1000 weakly 
interacting atoms) so that its three spin components, 
|0) = |F=l,mi. = 0) and |±1) = \F = l,mF = ±l), 
share the same spatial mode. In this so-called single- 
mode approximation (SMA) [2^, we can describe the 
spinor condensate (subject to a quadratic Zeeman effect) 
with the Hamiltonian [lal 

Hs =h\[{N+ - N^f + (27Vo - 1){N+ + iV_) 

+ 2clclc+c- + 2c^cLcoCo] + hq{N+ + 7V_), (1) 

where Ci is the field operator annihilating a bosonic atom 
in component with Ni = c\ci the corresponding atom 
number operator. A is a coefficient related to the spin- 
dependent part of the two-body interaction, and finally 
q the quadratic Zeeman shift. 

The cavity is assumed to support a 7r-polarized single 
traveling mode with frequency Wc, and is driven by an 
external probe field with an amplitude rj and frequency 
Wp. As in Ref. [11], both Wc and ujp are assumed to be 
sufficiently red-detuned from the F = 1 ^ F' = 1 atomic 
transition frequency uja so that the excited states can be 
adiabatically eliminated. Under such a circumstance, 
our cavity -t- condensate system (excluding the reservoir 
consisting of the cavity vacuum modes) is described by 
a total Hamiltonian H = Hg + where Hm is given 
explicitly by 



Hr, 



-h5a)a + hi] [a) + a) + hUo{N - Na)a^a, (2) 



with a being the field operator for annihilating a cav- 
ity photon and 5 = ujp — lUc the detuning of the probe 
relative to the cavity mode frequency. In addition to 
the cavity photon energy (the first term) and the Hamil- 
tonian simulating the process of pumping the cavity 
mode by the classical external probe field (the second 
term), a new term (the last term) appears in Eq. 
which characterizes the atom-photon interaction with 
an effective strength Uq = gg/ i^p " '^a) with go be- 
ing the atom-cavity mode coupling coefficient. Sev- 
eral comments arc in order. First, the selection rule 
for dipolc transitions involving 7r-polarized photons only 
permits the transitions between \F = l,mF = +1 (^1)) 



and \F' = l,m'p = +1 (—1)), and consequently, the last 
term, in the limit of far-off-resonant atom-photon inter- 
action, is expected to be proportional to , which 
is equivalent to iV — iVo when the definition for the total 
atom number, N = iV+ + iV_ -I- iVo, is taken into consid- 
eration. Second, one can express and in terms 
of No along with two constants of motion under the total 
Hamiltonian H: N and M = iV+ — iV_ (magnetization), 
and as a result, we will focus, from now on, our attention 
on the dynamics of (No)- Finally, we emphasize that 
the dispersive interaction term KUoNoo) a in Eq. ^ can 
cause the probe field to experience a phase shift propor- 
tional to (A^o), which, in the language of measurement 
in quantum optics, constitutes the (matter wave) signal 
we aim to determine from the measurement of the probe 
field. 

^ 
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FIG. 1: (color online) Schematic diagram of homodyne detection. 

These discussions lend itself naturally to the final com- 
ponent of our model. To begin with, we note that No does 
not commute with H and cannot serve as a quantum non- 



demolition measurement (QND) variable [21|, and thus 



the probe field is to remain as weak as the measurement 
permits in order to minimize its back-action on the sig- 
nal. As the photons leaking from the cavity are combined 
with those from a strong local oscillator prior to the mea- 
surement by the photodetector, the homodyne detection 
scheme illustrated in Fig.[T]can greatly enhance the signal 
to noise ratio while at the same time allows us to directly 
measure the quadrature phase amplitude and hence the 
(Ao)-depcndent probe phase shift as discussed above. In 
practice, in order to gain the spin dynamics, one must 
monitor the phase shift continually, and perform many 
runs of experiments, each of which provides a continual 
stream of information about (A^o)? before one can aver- 
age over all the runs to construct the ensemble average 
{{No)). The stochastic master equation approach, which 
combines the system-reservoir theory with the photon 
counting theory [l^], is believed to be an excellent tool 
to simulate such experimental processes. This is the ap- 
proach we take in the present study. 

We begin with the measurement outcome, i.e., pho- 
todetector current / (for a single run), which, after sub- 
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tracting the constant part due solely to the coherent local 
oscillator, reads [l^, [3 



2k (a + «■'■) + \/2k 



dW (t) 
dt ' 



(3) 



where k is the cavity decay rate, and dW{t)/dt represents 
Gaussian white noise with dW(t) an infinitesimal Wicn- 
ner increment satisfying the Ito rules, {{dW {t))) = and 
dW(tf = Hi. 

For the system subject to a continuous homodyne de- 
tection, its time evolution, conditioned on a given set 
of measurement outcomes, is described by the stochastic 
master equation (SME) 



dpc i.-fj 



2KC[d]pc 



'2k — - — H.[a\pc 



dt 



(4) 



where pc is the conditional density matrix operator for 
the cavity mode + condensate system, and C and % are 
the superoperators defined as 



L[x\p = xpx^ 



-x^xp 



-px^x, 



'H[x]p ~ xp + px^ — Tr (xp 



px'')p. 



The first term on the right-hand side of Eq. ^ represents 
the unitary evolution of the system under H. The sec- 
ond term describes the decay of the cavity, originating 
from coarse graining over the reservoir degrees of free- 
dom. The last term is related to the quantum state col- 
lapse accompanied by the detection of each photoelectron 
at the detector; the fact that it shares with the current 
in Eq. ([3|) the same noise term, dW{t)/dt, indicates that 
the evolution of pc is indeed conditioned on the current 
measurement. Both the second and the last term can 
affect the dynamics of the cavity field and the spinor 
condensate. 

To clearly show the measurement back-action on the 
spinor BEG, we consider that the measurement system 
operates in the regime where the cavity field decays at a 
rate much faster than the mean-field phase shift due both 
to the dispersive atom-photon coupling, and to the two- 
body s-wave scattering of atoms. Under such a condition, 
we can approximate a around a mean value a = {a} ~ 
■q/n (the field amplitude of an empty cavity) with a small 



fluctuation a'according to a = a 



and eliminate the 



modes defined by the bosonic operator a' adiabatically 
[igj . In this way, we arrive at the dimensionless SME for 
the conditional density operator psc = TrcavityPc of the 
spinor BEG alone 

^ = -i[H',Psc] + 2eC[No]psc + V2^^H[No]psc, 
dr dr 

(5) 

as well as the scaled photoelectric current 

/' = 2x/2e(iVo) + ^, (6) 

where t = \X\t is the scaled time, dW' (t) / dr the scaled 
white noise, and ^ = Ugr]/ \/ | A| the measurement 



strength [26|. If we were to ignore the measurement 
back-action, the spinor BEG would undergo a unitary 
evolution under the scaled effective Hamiltonian H' = 
[H, + hUn{N - No) \a\^ - M|a|^]/ft|A|. In what fol- 
lows, in order to highlight the essential physics, we fix 
the detuning to 5 = UqN so that H' becomes 

A 



+ 2cjcjc+c_ 4- 2c^c^CoCo] - q'Nq, 



(7) 



(after removing a constant term hqN) wnere q' = 
((7 + JJojap) / |A| is defined as a new quadratic Zeeman 
shift. 

The last two terms at the right-hand side of Eq. ([S]) 
represent the measurement back-action to the spinor con- 
densate. The first of these is proportional to the double 
commutator [TVq, [iVo,/3sc]], which represents a source of 
decoherence in the quantum dynamics. It tends to damp 
the off-diagonal elements of the density matrix under the 
basis of the measured observable Nq. It represents one 
form of measurement back-action as it originates from the 
fact that any measurements on the cavity mode require 
the use of an output coupler to couple the cavity mode 
to the field modes outside the cavity. The last term in 
Eq. (O can again be traced to the measurement induced 
state collapse in quantum mechanics, which is a stochas- 
tic process and hence depends on the white noise. The 
dynamics obtained directly from Eq. ([5]) is called the con- 
ditional dynamics, while that obtained after the ensem- 
ble average is called the deterministic dynamics. The last 
term of Eq. ^ therefore affects the conditional, but not 
the deterministic dynamics. Finally, in principle, there 
exists another type of measurement back-action - atom 
heating due to the fluctuation of the optical dipole force 
[2^. However, since such a fluctuation is proportional 
to the gradient of the cavity field intensity, we anticipate 
the heating effect in a traveling-wave cavity to be much 
weaker than what was observed in a standing-wave cav- 
ity [25} . and therefore we neglect it entirely in our work 
here. 



III. SPIN-MIXING DYNAMICS UNDER THE 
CONTINUOUS MEASUREMENTS 

A. Two-Atom Case 

In this section, we apply the formalism outlined in the 
previous section to a two-atom "toy" model, which, in 
principle, can be realized in optical lattices [l^, to il- 
lustrate the influence of measurement backaction. Due 
to the conservation of atom number and magnetization, 
a spin-1 BEG with two atoms and zero magnetization 
is effectively a spin- 1/2 system with two basis states 
|1) = |0,2,0) and |2) = |1,0, 1), where \N+,No,N^) is a 
Fock state with A^^ number of atoms in spin-i component. 
In this basis, H' (with q' = 0) has the following matrix 
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representation 



H' 



_ 2^/2 
\\\ V 2V2 -2 



(8) 



which has two eigenvalues, Ea = 2A/ |A| and Ef, = 
— 4A/|A|, and two corresponding eigenstates \a) = 
y273|l) + Vl73|2> and \b) = -^17311) + y^\2). 
Here, \a) is the ground state when A < (ferromag- 
netic case) and \b) is the ground state when A > 
(antiferromagnetic case). 

The dynamics of any atomic observable A in a partic- 
ular realization can be constructed with the help of SME 
dU) starting from 



(9) 



For our case here, we find that the dynamical equations 
for the three Hermitian operators defined as 



V 



(10) 







are closed and can be cast into a matrix form 
-4^2 
±2 -4^2 -p4^ 
±4^2 





2\/2e 



dW 
dr 



{Sy){S, 



(11) 



1 



Here the upper (lower) signs are for the antiferromag- 
netic (ferromagnetic) case. The terms associated with 
the coefficient — 4^^ represent the dampings, typical of 
the dynamics of an open system, which destroy the coher- 
ence and leaves the system in a mixed state composed of 
eigenstates of measured observables. In the current sec- 
tion, we only consider the antiferromagnetic case (A > 0) 
as it does not exhibit a qualitatively different dynamics 
from the ferromagnetic case. 

Consider two atoms that are initially prepared in state 
|1). Figures [2][a) and (b) illustrate, respectively, condi- 
tional and deterministic dynamics for systems with ^ = 
(red dotted lines), 0.5 (blue dashed lines), and 2 (black 
solid lines) . For the case of no measurement (^ = 0) , the 
population dynamics, {Nq)/N = {l + {Sz))/N, undergoes 
a Rabi-type oscillation with frequency \Ea — Eb\ = 6. For 
a relatively weak measurement (^ — 0.5), besides some 
superimposed noises, the oscillation in a single run be- 
gins to experience a diffusive phase shift relative to the 
one without measurement; this results in a damped oscil- 
lation that one expects when many oscillations with dif- 
ferent phase shifts are averaged. For a relatively strong 
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FIG. 2: (color online) The evolution of the fractional population in 
spin-0: Nq/N. (a) and (b) show the conditional and deterministic 
evolution, respectively, with the initial state |1) and q' = 0. The 
red dotted line is for non- measurement case ^ = 0, while the blue 
dashed hne and black solid hnc is for weak measurement = 0.5) 
and strong measurement = 2) case, respectively. In (c), the 
initial state is \b) (see text) and § = 0.1. The blue dashed line 
shows the conditional evolution while the black solid line shows 
the deterministic one. The red dotted line is for non-measurement 



measurement (^ = 2), the spin dynamics is drastically 
different. This is due to that the system is "watched" 
so frequently tha t Q uantum Zeno effect (QZE) begins 
to manifest itself [18[. Indeed, the conditional evolution 
indicates that repeated observations tend to trap the sys- 
tem in states |1) and |2), the only two states at which the 
noise term in Eq. (jlip vanishes. The extra time that the 
system spends either in |1) or |2) in the conditional evolu- 
tion slows down the transition of the initial state to other 
states as shown in the deterministic evolution, which de- 
cays exponentially without any oscillations. In both weak 
and strong measurements, deterministic evolutions con- 
verge to a mixed state {{Sx) = 0, {Sy) = 0, {Sz) = 0), the 
only fixed-point of the deterministic part of Eqs. (|lip 
at which the density matrix takes the diagonal form: 
p,c = 0.5|l)(l|+0.5|2) (2|. 

Let us now discuss Fig. [5] (c), which displays the dy- 
namics of two atoms initially prepared in the antiferro- 
magnetic ground state \b). In the absence of any measure- 
ments, as expected, the system stays in its ground state 
(red dotted line). For a weak measurement (^ = 0.1), 
the system develops a Rabi-type oscillations in the con- 



5 



ditional evolution, and is shown to attempt to con- 
verge to the mixed state in the deterministic evolution 
(black solid line). As before, QZE appears (not shown) 
when the measurement is sufficiently strong. In the two- 
atom case, a system starting from the antifcrromagnetic 
ground state exhibits similar dynamics as that starting 
from the ferromagnetic ground state. However, in the 
A^-atom case, as we show in the subsection below, due to 
the difference in the energy level structure and quantum 
fluctuation of TVq, the measurement back-action will have 
quite distinct effects on the ferromagnetic and antifcrro- 
magnetic ground states. 



B. A'^-Atom Case: Conditional Population 
Dynamics 

Now we illustrate the measurement back-action effect 
for a condensate with A = 100 atoms. We also adopt re- 
alistic parameters: k = 27rx 100 MHz, go — 27rx 1.6 MHz, 
and A = 27r X 20 Hz for sodium atoms with a typical den- 
sity lO^"' cm^'^ [lo|. We estimate that the value of ^ lies 
in the range between 10~^ and 10~^ which means mea- 
surements here are always very weak. It is impossible to 
find a set of observables that are closed under Eq. ^ as 
in the two-atom case. However, under the assumption of 
perfect detection (with unit detection efficiency), we can 
unravel SME ([5]) into a equivalent stochastic Schrodinger 
equation (SSE) [g^] in the sense psc = K'sc) (V'sd as 



dW 



-V2e(iVo- (TVo))- 



dt 



(12) 



which allows the dynamics of any observables to be cal- 
culated exactly. The simulation is performed using a 
fourth-order Runge-Kutta method for the deterministic 
part, and a first-order stochastic Runge-Kutta method 
for the noise part. Furthermore, the SSE shows more 
clearly that the measurement back-action effects are de- 
pendent not only on the measurement strength, but also 
on the quantum fluctuation of the measured observable 
Ao. 

For the antifcrromagnetic case, the ground state for 
H' (with q' = 0) is unique and is given by a superpo- 
sition of all the Fock states in which the spin-1 and —1 
components share the same atom number: 

lN/2] 

= J2 A,\k,N-2k,k), 

k=0 

where the amplitudes Ak obey the recursion relation Q 





2/cH 


-2 




2k -\ 


- \ 



A 



k~l- 



and the number-fluctuation in spin-0 component is super- 
Poissonian for A^ ^ 1 This indicates that the anti- 
ferromagnetic ground state will be quite sensitive to the 
measurement back-action effect. 
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In this state, the average atom numbers in each spin com- 
ponent are equal, i.e., {N+} = (A_) = (A"o) = A^/3, 
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FIG. 3: Left column: Evolution of normalized particle number in 
spin-0 component under a single run of measurements with different 
strength (a) ^ = 0.001, (b) 0.01, and (c) 0.1. The initial state is 
the anti-ferromagnetic ground state (see the text) with total atom 
number N = 100. The corresponding Fourier spectra are shown in 
the right column. 

Figures |3] illustrate the conditional population dy- 
namics of atoms initially prepared in the antifcrromag- 
netic ground state and the corresponding Fourier spectra 
for various measurement strengths. As Fig. EJa) illus- 
trates, a measurement as weak as ^ = 0.001 can induce 
the system to oscillate predominantly at frequency 6A 
(Fig. |31[d)), the ffist excited frequency of the many-body 
system described H' in Eq. ([7]). As the measurement 
strength increases, more and more high-frequency com- 
ponents contribute to the evolutions. Fig. ^h) is pro- 
duced with ^ = 0.01 which is ten times stronger than in 
Fig. [31Ja). Indeed, instead of one peak, its Fourier spec- 
trum [Fig. ^e)] displays four peaks corresponding to the 
first fourth excited frequencies. Increasing ^ by another 
factor of ten to ^ = 0.1 leads to a more chaotic evolution 
as con&mcd both by the population dynamics in Fig. 
ISjc) and by the corresponding Fourier spectrum in Figl3] 
(f). Here, as a result of a dramatic increase in the num- 
ber of eigenstates to which the system can collapse, QZE 
becomes more complicated. In Fig. ^c), only the tran- 
sitions to \N/2,0,N/2) and |0,Ar,0) are (dimly) visible 
because these two have the smallest transition moments 
to their neighboring eigenstates. In order to see QZE 
involving other transitions, we find from our numerical 
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simulations that strong measurements with ^ > 10 are 
typically needed. 




X (n/|A,| 



FIG. 4: Same as Fig. [S] except that for a ferromagnetic conden- 
sate. Insets show details of the oscillations. 



The measurement back-action on ferromagnetic spinor 
condensates (A < 0) will take somewhat different effects 
because their distinct energy-level structures and quan- 
tum statistical properties. The ferromagnetic ground 
state is 27V + 1-fold degenerate, which reads Q 

IV-m) = ^fc \k,N-2k-m,k + m), 

k 

where m = 0, ±1, . . . , ±iV. Contrary to the anti- 
ferromagnetic ground state, these states possess sub- 
Poissonian fluctuations in Nq. To demonstrate the back- 
action, we choose the one with m = as the initial state, 
which has the largest fluctuation in Nq. 

In Fig. |4l^a), the atomic population exhibits a weak 
Rabi-type oscillations with an amplitude much smaller 
than that for the antiferromagnetic case under the same 
measurement strength. The main reason for this reduc- 
tion is, as indicated by the Fourier spectrum in Fig. 
IDJd), that the first excited frequency is located around 
cj = 398, which is much higher than that for the ferro- 
magnetic case and hence is much difficult to excite. The 
reduction in the variance of A^o may also weaken the ef- 
fect of measurement back-action. With the increase of 
the measurement strength, similar to the antiferromag- 
netic case, the spin populations oscillate with multiple 
frequencies and become irregular with some evidence of 
QZE, as shown in Figs. |li;b)-(f). 



C. A'^-Atom Case: Comparison between 
Conditional and Deterministic Population Dynamics 

The two examples considered in the previous subsec- 
tion demonstrate the effect of the measurement back- 
action on the population dynamics of a single experimen- 
tal realization, and as in the two-atom case, the determin- 
istic dynamics will emerge from the average over many 
runs of numerical simulations. In order to make a smooth 
transition to the subject discussed in the next section, in- 
stead of pursuing such simulations with the two examples 
considered above, we seek to show the effect of averaging 
over many runs from a spinor BEC initially prepared in 
its mean- field ground state ji/; (r = 0)) = |0, TV, 0), where 
all the atoms reside in the spin-0 component. 
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FIG. 5: (color onhne) (a) Time evolutions of (No)/N for q' = 10 
(upper curves) and q' = (lower curves) case with the same ini- 
tial state |0, A'^, 0). The blue curves correspond to the evolution 
without measurement. The gray dotted curves display a variety 
of evolutions conditioned on measurement outcomes with measure- 
ment strength ^ = 0.01, and the red curves are given by averaging 
over 10 conditional evolutions, (b) The number fluctuation of spin- 

component Var(7Vo)=-y/ (Nq) — (Nq)^ , with solid and dashed line 
corresponding to q' = 10 and q' = case, respectively. 

The spin-mixing dynamics, in the absence of the probe, 
can be well understood from the quantum-fluctuation- 
driven harmonic oscillator model [15| . In Fig. [SJa), we 
compare the spin dynamics between q' ~ 10 (upper solid 
blue curve) and q' — (low solid blue curve). In the 
former case, the oscillations are weak and approximately 
harmonic, while the latter case exhibits oscillations that 
are clearly of anharmonic nature. This is because quan- 
tum fluctuation in Nq is larger for small q' than for large 
q' as shown in Fig. EJb). The spin dynamics for q' ~ 
in a longer time scale is shown by the dashed blue curve 
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in Fig. [6ja) , which clearly demonstrates a typical quan- 
tum behavior - collapse and revival of spin oscillations. 
The particle number distribution quickly collapses to a 
metastable regime with (iV+) = {N_) = {Nq)/2 = N/i 
after a time Tc — {2^/N)~^. This metastable regime is 
followed by several spin oscillations and the cycle repeats 
itself at a time interval = tt Q . 




0.0-^ ^ , ^ , ■ , ■ , ■ 1 X 

2 4 6 8 10 



FIG. 6: (color online) Comparing between the evolution without 
measurement (a) and the deterministic evolution (b) for the q' = 
case in a long time scale. 

In the presence of the probe, the spin-mixing dynamics 
will be affected by the measurement back-action. Var- 
ious conditional evolutions with an intermediate mea- 
surement strength ^ = 0.01 are plotted in gray dotted 
curves in Fig. [Sl^a). For the case q' = 10, they are not 
much different from the non-measurement evolution, so 
that averaging over 10 conditional evolutions appears suf- 
ficient to reveal the deterministic spin dynamics. The 
quantum measurement back-action are restrained by the 
small quantum fluctuation. In contrast, for the case 
q' ~ 0, they differ from the non-measurement evolu- 
tion quite appreciably. In this case, averaging over 10 
conditional evolutions (solid red curve) cannot produce 
the anticipated deterministic dynamics and the match is 
particularly poor in the metastable regime due to the 
large quantum fluctuation there. It requires more runs 
of measurements to reveal the deterministic spin evolu- 
tion. The curve in Fig. IHl^b) represents its deterministic 
evolution given by averaging over 100 conditional evo- 
lutions, which in a short time scale traces out the an- 
harmonic spin oscillations clearly but indicates that the 
oscillations are gradually damped for a long time evo- 
lution. The damping rate is proportional to the mea- 
surement strength. Finally, the BEC converges to a 
mixed state characterized with a diagonal density matrix, 

Psc = J2'i!ilPk\k,N-2k,k){k,N-2k,k\, where the 
probability distribution function Pk is found (not shown) 
to be a constant independent of fc or = 1/ {N/2 -I- 1) to 
be precise. All these are due to the decoherence induced 



by measurement as discussed in the two-atom case. 



IV. THE MEASUREMENT OUTCOME: 
PHOTOELECTRIC CURRENT 

The numerical simulations we have considered so far 
show that although each run results in a different condi- 
tional evolution (Aq), an ensemble average over dozens 
of these runs can already capture quite well the deter- 
ministic quantum spin-mixing dynamics. However, what 
is accessible in experiments is not (A^o) but the photo- 
electric current [Eq. Thus, in practice, {{No)) must 
be inferred by averaging the current over many runs of 
measurements. As it turns out, it requires far more runs 
to reveal ((Aq)) indirectly from the ensemble average of 
the photoelectric current than directly from the ensemble 
average of conditional population dynamics. 

Figures [7] (a) and (b) show the ensemble averages of 
the photodetector current /' for A^ = 10 and A^ = 100, 
respectively. The black curves represent the results given 
by 100 runs of measurements. As can be seen, it is vir- 
tually impossible to extract the deterministic evolutions 
of {{Nq)) (shown in the insets) as they are dominated 
by white noise. In principle, one can suppress the noise 
by averaging over more and more currents; this is evi- 
dent from the examples obtained when we increase the 
number of runs to 10'^ (blue curve) and then to 10^ (red 
curve). However, only in A^ = 100 case can the spin 
oscillations of deterministic nature be (vaguely) recog- 
nized. As for A^ = 10 case, averaging the current over 
10^ runs of measurements still do not allow us to extract 
the signal. 

It appears that one could increase the measurement 
strength ^ instead of the number of runs to enhance the 
signal to noise ratio according to Eq. ([5]). However, in 
quantum measurements, the system dynamics is condi- 
tioned on the detection outcomes; increasing the mea- 
surement strength also enhances the quantum measure- 
ment back-action. First, according to the discussion in 
Sec. II, strong measurement renders a large decoherence 
to the measured quantum state so that the spin oscilla- 
tions is rapidly damped. Second, an increase in ^ will 
increase the white noise in the stochastic Schrodinger 
equation, which in turn demands more runs to recover 
the deterministic dynamics. Thus, there is a limitation 
to what we can do to improve the signal to noise ratio by 
increasing the measurement strength. 

An alternative is to increase the atom number. Not 
only does it enhance the signal part of the current in 
Eq. but also it reduces the quantum fluctuation of 
A^o and thus the related measurement back-action. The 
net effect is the reduction in the number of required runs. 
But, as the atom number increases, the mean-field dy- 
namics will gradually dominate [1], defeating the goal 
of extracting beyond-mean-ficld quantum dynamics from 
this measurement scheme. A possible way to increase the 
signal to noise ratio without raising the atom number is 
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the stochastic master equation approach, we have per- 
formed a detailed study of the quantum measurement 
back-action on the spin population dynamics in the bad 
cavity limit. We have used a simple two-atom system to 
illustrate both the measurement-induced quantum Zeno 
effect and the measurement-induced diffusive quantum 
dynamics. We have applied the physical intuitions gained 
from the two-atom model to understand the measure- 
ment back-action on the spin population dynamics in a 
spin-1 BEC. We have shown that the effect of back-action 
is sensitive to the quantum fluctuation of the spinor con- 
densate. Finally, we stress that this study is motivated 
by recent proposals for using measurement techniques 
popular in cavity quantum optics to probe the quantum 
dynamics of small condensates. An important point we 
aim to make in this work is that when applying opti- 
cal detection techniques to small condensates, one needs 
to pay close attention to quantum fluctuations, which 
are typically ignored for large condensates. Indeed, we 
have shown that due to the quantum measurement back- 
action, the number of runs of measurements, needed to 
recover the deterministic population dynamics from the 
ensemble average of the photoelectric current, increases 
as the number of atoms in the condensate decreases, sug- 
gesting that the scheme is not practical for sufficiently 
small condensates where the number of runs can become 
unrealistically large. 



FIG. 7: (color online) Measurement outeomes given by averaging 
over 10^ (black), 10"^ (blue), and 10^ (red) photodetector currents. 
The total atom number TV is 10 in (a) and 100 in (b), and other 
parameters are the same as in Fig. [5] The insets show the deter- 
ministic evolution of (No)- 

to process the current signal using methods such as filter- 
ing high frequency components and averaging over sliding 
windows [28|. However, our study shows that more than 
10'^ runs are still needed before we can reveal the spin 
dynamics for a small BEC with less than 100 atoms. 

V. CONCLUSION AND REMARKS 

In this work, we have considered a homodyne detec- 
tion scheme, which is designed to make a continuous 
measurement of the quantum spin-mixing dynamics of 
a small F=l spinor BEC inside an optical cavity. Using 
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